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ABSTRACT 

Magnetars are neutron stars possessing magnetic field of about lO^'^'^-lO"'^^ G at the surface. Ther- 
modynamic properties of neutron star matter approximated by pure neutron matter are considered at 
finite temperature in strong magnetic fields up to 10^* G which could be relevant for the interior regions 
of magnetars. In the model with the Skyrme effective interaction, it is shown that a thermodynami- 
cally stable branch of solutions for the spin polarization parameter corresponds to the case when the 
majority of neutron spins are oriented opposite to the direction of the magnetic field (negative spin 
polarization). Besides, the self-consistent equations, beginning from some threshold density, have also 
two other branches of solutions corresponding to positive spin polarization. The influence of finite tem- 
peratures on spin polarization remains moderate in the Skyrme model up to temperatures relevant for 
protoneutron stars. In particular, the scenario with the metastable state characterized by positive spin 
polarization, considered at zero temperature in Phys. Rev. C 80, 065801 (2009), is preserved at finite 
temperatures as well. It is shown that above certain density the entropy for various branches of spin 
polarization in neutron matter with the Skyrme interaction in a strong magnetic field demonstrates 
the unusual behavior being larger than that of the nonpolarized state. By providing the corresponding 
low-temperature analysis, it is clarified that this unexpected behavior should be addressed to the de- 
pendence of the entropy of a spin polarized state on the effective masses of neutrons with spin up and 
spin down, and to a certain constraint on them which is violated in the respective density range. 

Key words : Neutron star models, magnetar, neutron matter, Skyrme interaction, strong magnetic 
field, spin polarization, finite temperature 



I. INTRODUCTION 

Magnetars are strongly magnetized neutron stars 
(Duncan & Thompson 1992) with the magnetic field 
strength at the surface of about 10^'*-10^^ G (Thomp- 
son & Duncan 1996, Ibrahim et al. 2002). Such huge 
magnetic fields can be inferred from observations of 
magnetar periods and spin-down rates, or from hydro- 
gen spectral lines. Among possible classes of various 
neutron stars, soft gamma-ray repeaters and anoma- 
lous X-ray pulsars are believed to be most probable 
candidates for these ultrastrong magnetized astrophys- 
ical bodies (Woods & Thompson 2006). Magnetars are 
relatively frequent objects in the Universe and comprise 
about 10% of the whole population of neutron stars 
(Kouveliotou et al. 1998). In the interior of a magnetar 
the magnetic field strength may be even larger, reach- 
ing values of about 10^^ G (Chakrabarty, Bandyopad- 
hyay & Pal 1997, Broderick, Prakash & Lattimer 2000). 
Therefore, magnetars provide a unique playground for 
studying properties of neutron star matter under ex- 
treme conditions of density and magnetic field strength 
(Chakrabarty, Bandyopadhyay & Pal 1997, Broderick, 
Prakash & Lattimer 2000, Cardall, Prakash & Lattimer 



2001, Perez-Garcia 2008, Isayev & Yang 2009), which 
are inaccessible in the terrestrial laboratories. 

In the recent study by Perez-Garcia (2008), neutron 
star matter was approximated by pure neutron mat- 
ter in a model with the effective nuclear forces. It has 
been shown that the behavior of spin polarization of 
neutron matter in the high density region in a strong 
magnetic field crucially depends on whether neutron 
matter develops a spontaneous spin polarization (in 
the absence of a magnetic field) at several times nu- 
clear matter saturation density, or the appearance of a 
spontaneous polarization is not allowed at the relevant 
densities (or delayed to much higher densities). The 
first case is usual for the Skyrme forces (Rice 1969, Sil- 
verstein 1969, 0stgaard 1970, Viduarre, Navarro & 
Bernabeu 1984, Reddy et al. 1999, Akhiezer, Laskin 
& Peletminsky 1996, Marcos et al. 1991, Maruyama 
& Tatsumi 2001, Beraudo et al. 2004, Kutschera 
& Wojcik 1994, Isayev 2003, Isayev & Yang 2004a, 
Rios, Polls & Vidafia 2005, Isayev 2006), while the 
second one is characteristic for the realistic nucleon- 
nucleon (NN) interaction (Pandharipande, Garde & 
Srivastava 1972, Backmann & Kallman 1973, Haensel 
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1975, Vidana, Polls & Ramos 2002, Fantoni, Sarsa, 
& Schmidt 2001, Sammarruca & Krastev 2007, Bor- 
dbar & Bigdeli 2007). In the former case, a ferro- 
magnetic transition to a totally spin polarized state 
occurs while in the latter case a ferromagnetic transi- 
tion is excluded at all relevant densities and the spin 
polarization remains quite low even in the high den- 
sity region. If a spontaneous ferromagnetic transition 
is allowed, it was shown in the subsequent model con- 
sideration with the Skyrme effective forces ( Isayev & 
Yang 2009) that the self-consistent equations for the 
spin polarization parameter at nonzero magnetic field 
have not only solutions corresponding to negative spin 
polarization (with the majority of neutron spins ori- 
ented opposite to the direction of the magnetic field) 
but, because of the strong spin-dependent medium cor- 
relations in the high-density region, also the solutions 
with positive spin polarization. In the last case, the 
formation of a metastable state with the majority of 
neutron spins oriented along the magnetic field is pos- 
sible in the high-density interior of a neutron star. 

In this study, we extend our previous consideration 
of neutron matter with Skyrme forces in a strong mag- 
netic field, given at zero tcmpcratiirc (Isayev & Yang 
2009), to finite temperatures up to a few tens of MeV 
being relevant for protoneutron stars. As a frame- 
work for consideration, we use a Fermi liquid approach 
for the description of nuclear matter (Akhiezer et al. 
1994, Akhiezer et al. 1997, Isayev & Yang 2006). One 
of the goals of the research is to study the impact of 
finite temperatures on spin polarized states in neutron 
matter in a strong magnetic field, and, in particular, 
to clarify whether the formation of a metastable state 
with positive spin polarization is possible in the high- 
density region of neutron matter. Besides, we provide 
a fully self-consistent finite temperature calculation of 
the basic thermodynamic functions of spin polarized 
neutron matter in a strong magnetic field. It will be 
shown that beginning from some density the entropy for 
various branches of spin polarization in neutron matter 
in a strong magnetic field demonstrates the unusual be- 
havior being larger than that of nonpolarized neutron 
matter. We relate this unexpected result to the depen- 
dence of the entropy on the effective masses of neutrons 
with spin up and spin down, and to the violation of a 
certain constraint on them in the corresponding density 
range. 

Note that we consider thermodynamic properties of 
spin polarized states in neutron matter in a strong mag- 
netic field up to the high density region relevant for as- 
trophysics. Nevertheless, we take into account the nu- 
cleon degrees of freedom only, although other degrees 
of freedom, such as pions, hyperons, kaons, or quarks 
could be important at such high densities. 

II. BASIC EQUATIONS 

Here we present only the basic formulae necessary for 
further calculations, although more details concerning 



a Fermi-liquid approach to neutron matter in a strong 
magnetic field can be found in our earlier work (Isayev 
& Yang 2009). The normal (nonsuperfiuid) states of 
neutron matter are described by the normal distribu- 
tion function of neutrons /k^kj — Tr f^a+a^i, where 
K = (p, (t), p is momentum, a is the projection of spin 
on the third axis, and g is the density matrix of the 
system. Further it will be assumed that the third axis 
is directed along the external magnetic field H. Given 
the possibility for alignment of neutron spins along or 
opposite to the magnetic field H, the normal distribu- 
tion function of neutrons and single particle energy can 
be expanded in the Pauli matrices crj in spin space 

/(p) = /o(p)ao + /3(p)a3, (1) 

£(p) = eo(p)o-o + e3(p)o-3- 

Expressions for the distribution functions /o,/3 in 
terms of the quantities s read (Isayev 2003, Isayev & 
Yang 2004a) 

/o = i{n(u;+)+n(u;_)}, (2) 
/3 = \{n{^+)-n{uJ-)}. 

Here n(w) = {exp(Fo'^) + and 

w± = Co ±6, (3) 
6 = £0 - Ato, 6 = ^3, 

Ho being the chemical potential of neutrons. The 
branches (jj± of the quasiparticle spectrum correspond 
to neutrons with spin up and spin down. 

The distribution functions / should satisfy the norma- 
lization conditions 

^E/o(P) = (4) 
p 

2 X - 

yZ^faip) = Qt-Qi = ^g- (5) 
p 

Here Q = + Qi is the total density of neutron mat- 
ter, and g^ are the neutron number densities with 
spin up and spin down, respectively. The quantity Ag 
may be regarded as the neutron spin order parame- 
ter. It determines the magnetization of the system 
M = /i„A£i, being the neutron magnetic moment. 
The magnetization may contribute to the internal mag- 
netic field B = H + A-kM . However, as we discussed 
earlier (Isayev & Yang 2009), and, analogously to the 
previous works (Perez-Garcia 2008, Broderick, Prakash 
& Lattimer 2000), we will assume that the contribu- 
tion of the magnetization to the magnetic field B re- 
mains small for all relevant densities and magnetic field 
strengths, and, hence, B k, H. 
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The self-consistent equations for the components of 
the single-particle energy have the form (Isayev & Yang 
2009) 

^o(p) = £o(p) +£'o(p) -Mo, (6) 
6(P) = -finH + ssip). (7) 

2 

Here eo(p) = is the free single particle spectrum, 
mo is the bare mass of a neutron, and £o , £3 are the 
Fermi liquid (FL) corrections to the free single par- 
ticle spectrum, related to the normal FL amplitudes 
[/o"(k),[/f (k) by formulas 



^~o(p) = ^E^o(k)/o(q), k: 



p q 



£3(p) = ^Ef^rw/3(q)- 



(8) 
(9) 



To obtain numerical results, we utilize the effective 
Skyrme interaction (Vautherin & Brink 1972). Expres- 
sions for the normal FL amplitudes in terms of the 
Skyrme force parameters were written by Akhiezer et 
al. 1997, Isayev & Yang 2006. Thus, using expres- 
sions (2) for the distribution functions /, we obtain the 
self-consistent equations (6), (7) for the components of 
the single- particle energy ^0 (p) and ^3 (p) , which should 
be solved jointly with the normalization conditions (4), 
(5). 

Note that spin ordering in neutron matter can be 
characterized by the neutron spin polarization param- 
eter 

^ gt - Ql ^ 

g ~ g ' 

The number densities of neutrons with spin up and spin 
down are related to the spin polarization parameter 11 
by formulas 

^>t = |(i + n), ^,; = |(i-n). (10) 

To examine the thermodynamic stability of different 
solutions of the self-consistent equations, it is necessary 
to compare the corresponding free energies F = E—TS, 
where the energy functional E is characterized by two 
FL amplitudes Ul^,U^ (Isayev & Yang 2006) and the 
entropy reads 



(11) 



p a= 



+n{u!a-) lnn{uja-)}, n{uj) = 1 — ■n{uj). 



III. 



SOLUTIONS OF SELF-CONSISTENT 
EQUATIONS AT FINITE T. THERMO- 
DYNAMIC FUNCTIONS 



The self-consistent equations were analyzed at zero 
temperature by Isayev & Yang 2009 for the magnetic 
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Fig. 1. — (Color online) Neutron spin polarization pa- 
rameter as a function of density at T = 20 MeV and dif- 
ferent magnetic field strengths for SLy7 interaction. The 
branches of spontaneous polarization II;^ , IIil are shown by 
solid curves. 



field strengths up to Hmax ^ 10^^ G, allowed by a 
scalar virial theorem (Lai & Shapiro 1991), in the 
model consideration with SLy4 and SLy7 Skyrme ef- 
fective forces (Chabanat et al. 1998). These Skyrme 
parametrizations were constrained originally to repro- 
duce the results of microscopic neutron matter calcu- 
lations (pressure-versus-density curve) and give neu- 
tron star models in a broad agreement with the ob- 
servables such as the minimum rotation period, grav- 
itational mass-radius relation, the binding energy, re- 
leased in supernova collapse, etc. (Rikovska Stone et al. 
2003). It was shown that a thermodynamically stable 
branch of solutions for the spin-polarization parameter 
as a function of density corresponds to negative spin 
polarization when the majority of neutron spins are 
oriented opposite to the direction of the magnetic field. 
Besides, beginning from some threshold density depen- 
dent on the magnetic field strength, the state with pos- 
itive spin polarization can be realized as a metastable 
state in the high-density region of neutron matter in 
the model with the Skyrme effective forces. 

Here we study the impact of finite temperatures on 
spin polarization of neutron matter in a strong mag- 
netic field. In particular, we are interested to learn 
whether the metastable state with positive spin polar- 
ization will survive at temperatures about several tens 
of MeV relevant for protoneutron stars. To that end, we 
directly find solutions of the self-consistent equations 
at nonzero temperature. Because the results of cal- 
culations with SLy4 and SLy7 Skyrme forces are very 
close, here we present the obtained dependences only 
for the SLy7 Skyrme interaction. 

Fig. 1 shows the neutron spin polarization parame- 
ter as a function of density, normalized to the nuclear 
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saturation density go (for SLy7 force, go — 0.158 fm"^), 
for a set of fixed values of the magnetic field strength at 
the temperature T = 20 MeV. The branches of sponta- 
neous polarization n^,n(^, corresponding to the van- 
ishing magnetic field, are shown by solid curves. As 
one can see, the obtained dependences qualitatively are 
similar to those obtained at zero temperature . The 
branches of spontaneous polarization are modified dif- 
ferently by the magnetic field: The branch 11,^ (gi) turns 
to the branch Ili{g) with negative spin polarization 
while the branch 11 J (p) splits into two branches, I12{q) 
and Ils^g), corresponding to positive spin polarization. 
For the lower branch Ili{g), there are three charac- 
teristic density domains. At low densities g ^ 0.5 go, 
the magnitude of the spin polarization parameter in- 
creases with decreasing density. At intermediate den- 
sities 0.5 go < g <^ 3 go, there is a plateau in the Tli(g) 
dependence, whose characteristic value depends on H 
(at the given temperature). At densities g ^ 3 go, the 
absolute value of the spin polarization parameter in- 
creases with density and tends to unity. Note that in 
the low-density domain the possibility of the appear- 
ance of a "nuclear magnetic pasta" and its impact on 
the neutrino opacities in the protoneutron star early 
cooling stage should be explored in a more elaborated 
analysis as discussed in detail by Perez-Garcia (2008). 

The upper branches 112 (p) and H^lg), correspond- 
ing to positive spin polarization, appear stepwise at 
the same threshold density gth dependent on the mag- 
netic field (at the given temperature) and being larger 
than the critical density of spontaneous spin instability 
in neutron matter. For the branch Il2{g), the spin po- 
larization parameter decreases with density and tends 
to zero while for the branch Hsig) it increases with 
density and approaches unity. Because of the nega- 
tive value of the neutron magnetic moment, the mag- 
netic field tends to orient the neutron spins opposite 
to the magnetic field direction. As a result, the spin 
polarization parameter for the branches Il2{g), Tlsig) 
with positive spin polarization is smaller than that for 
the branch of spontaneous polarization n,^, and, vice 
versa, the magnitude of the spin polarization parame- 
ter for the branch Ili{g) with negative spin polarization 
is larger than the corresponding value for the branch of 
spontaneous polarization IIq . 

Thus, at densities larger than gth, we have three 
branches of solutions: one of them, Ili{g), with neg- 
ative spin polarization and two others, 112 (j?) and 
H^lg), with positive polarization. In order to clar- 
ify which branch is thermodynamically preferable we 
should compare the corresponding free energies. Fig. 2 
shows the free energy per neutron as a function of 
density at T = 20 MeV and H = lO^^ G for these 
three branches, compared with the free energy per neu- 
tron for a spontaneously polarized state [the branches 
n^(p)]. It is seen that the state with the majority 
of neutron spins oriented opposite to the direction of 
the magnetic field [the branch ni(£<)] has a lowest free 
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Fig. 2. — (Color online) Free energy per neutron as a 
function of density at T = 20 MeV for different branches 
X\\{^q)-X\z(,q) of solutions of the self-consistent equations at 
H = 10^* G for SLy7 interaction, including a spontaneously 
polarized state. 



energy. This result is intuitively clear, since magnetic 
field tends to direct the neutron spins opposite to H, as 
mentioned earlier. However, the state, described by the 
branch ^-ii^g) with positive spin polarization, has the 
free energy very close to that of the thermodynamically 
stable state. This means that despite the presence of a 
strong magnetic field H ^ 10^^ G, the state with the 
majority of neutron spins directed along the magnetic 
field can be realized as a metastable state in the dense 
core of a neutron star in the model consideration with 
the Skyrme effective interaction. In this scenario, be- 
cause such states exist only at densities g > gth, under 
decreasing density (going from the interior to the outer 
regions of a magnetar) a metastable state with positive 
spin polarization at the threshold density gth changes 
to a thermodynamically stable state with negative spin 
polarization. 

Note at this point that the possibility of the appear- 
ance of a metastable state with the majority of neu- 
tron spins oriented along the magnetic field was missed 
in the finite temperature calculations by Perez-Garcia 
(2008) because only one branch of solutions for the spin 
polarization parameter, corresponding to the thermo- 
dynamically stable state, was found there for neutron 
matter in a strong magnetic field with the same Skyrme 
interaction. 

An unexpected moment appears if we consider sep- 
arately the entropy for various branches of spin polar- 
ization. Fig. 3 shows the difference between the en- 
tropy per neutron for the branches Hi -Ha and that of 
the nonpolarized state (with 11 = at iJ = 0) as a 
function of density at the fixed magnetic field strength 
H = 10^^ G and different temperatures. Contrarily to 
the intuitively expected behavior, the entropy for the 
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Fig. 3. — (Color online) The entropy per neutron mea- 
sured from its value in the nonpolarized state for the IIi- 
Ds branches of spin polarization as a function of density 
at = 10^* G and different temperatures for the SLy7 
interaction. 



branch Hi, beginning from a certain density weakly 
dependent on temperature (at the given magnetic field 
strength), is larger than the entropy of the nonpolar- 
ized state. Besides, the entropy for the branches 112 
and 113 is larger than that for the nonpolarized state 
for all relevant densities. It looks like the spin polarized 
states described by the Hi -Ha branches of spin polar- 
ization are less ordered than the nonpolarized state for 
the corresponding densities. 

In order to understand qualitatively such an ex- 
pected behavior, let us consider the low-temperature 
expansion for the entropy in terms of the effective 
masses of spin- up and spin-down neutrons. The density 
of entropy (11) can be written in the form 



1 










H 






13 



E 



(12) 



T 



where 



J Art) 



-dx 



is Fermi-Dirac integral of the order [la is the ef- 
fective chemical potential of neutrons with spin up 
[a = V) and spin down(cr = — ), whose explicit ex- 
pression was written by Isayev & Yang (2009). In 
the low-temperature limit, 770- ^ 1, and, providing 
the corresponding expansion of Fermi-Dirac integrals 



in Eq. (12), one gets 



where e„ ~ ^ 



2£, 



-T, 



(13) 



— is the Fermi energy of neutrons with 

spin up and spin down, ka = {^-k^ QaY^^ being the re- 
spective Fermi momentum. The low-temperature ex- 
pansion (13) is valid till T/epa ^ 1. Then, requiring 
for the difference between the entropies of spin polar- 
ized and nonpolarized states to be negative, one gets 
the constraint on the effective masses Wnt ^^.d m„j^ of 
neutrons with spin up and spin down in a spin polarized 
state: 

i5 = !!^(i + n)^ + !^(i-n)^-2<o, (u) 

m„ nrin 

Here m„ is the effective mass of a neutron in nonpolar- 
ized neutron matter (Isayev & Yang 2004a), 



2rn„ 2too 



+ f[ti(l-a;i) + 3^2(1 + 2:2)], (15) 



and expressions for the effective masses m„-f- , m„4, read 
(Isayev & Yang 2009) 



2m 



til) 



In Eqs. (15), (16), ti and Xj are some phenomenolog- 
ical parameters, specifying a given parametrization of 
the Skyrme interaction. According to Eq. (10), the 
number densities of neutrons with spin up and spin 
down are determined by the spin polarization parame- 
ter n, which, after the self-consistent determination, is 
a function of the thermodynamic parameters g, T, H. 
Hence, the condition (14) determines the region in the 
domain of admissible values of the thermodynamic pa- 
rameters, where the entropy of a spin polarized state 
demonstrates the expected behavior. 

Fig. 4 shows the left-hand side D of constraint (14) 
for the branches H1-H3 of spin polarization as a func- 
tion of density at H = 10^^ G and temperature T = 
5 MeV, at which the condition ep^/T ^ 1 holds true. 
For the branch Hi, the difference D is negative up to 
the density gs « 0.41go (see the insert in Fig. 4), where 
it changes sign and remains positive for all larger densi- 
ties. Besides, for the branches H2 and H3, the difference 
D is positive for all densities at which the corresponding 
solutions of the self-consistent equations exist. These 
features explain the above mentioned unusual behavior 
of the entropy for various branches of spin polarization 
in neutron matter with the Skyrme interaction under 
the presence of a strong magnetic field. Note that such 
an unexpected behavior of the entropy was also found 
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Fig. 4. — (Color online) The difference D in constraint (14) 
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interaction. 



for the states with spontaneous spin polarization (in the 
absence of the magnetic field) in neutron matter with 
the Skyrme interaction (Rios, Polls & Vidafia 2005) and 
antiferromagnetically spin ordered states in symmetric 
nuclear matter with the Gogny DIS interaction (Isayev 
& Yang 2004b, Isayev 2005, Isayev 2007). 

IV. CONCLUSIONS 

We have studied the impact of finite temperatures 
on the spin structure in the magnetar interior, approxi- 
mating the neutron star matter by pure neutron matter 
and taking the Skyrme effective interaction as a poten- 
tial of NN interaction (SLy7 parametrization) . Accord- 
ing to the scalar virial theorem, strong magnetic fields 
up to 10^® G can be relevant for the interior regions of 
magnetars. It has been shown that, together with the 
thermodynamically stable branch of solutions for the 
spin polarization parameter corresponding to the case 
when the majority of neutron spins are oriented op- 
posite to the direction of the magnetic field (negative 
spin polarization), the self-consistent equations, begin- 
ning from some threshold density, have also two other 
branches of solutions corresponding to positive spin po- 
larization. The influence of finite temperatures on spin 
polarization remains moderate in the Skyrme model, 
at least, up to temperatures relevant for protoneu- 
tron stars. In particular, a thermodynamic analysis, 
based on the calculation of the free energy for different 
branches of spin polarization, shows that the scenario 
with the metastable state characterized by positive spin 
polarization, considered at zero temperature (Isayev & 
Yang 2009), is preserved at finite temperatures as well. 
This is one of the important conclusions of the given 
research. The possible existence of a metastable state 



with positive spin polarization will affect the neutrino 
opacities of a neutron star matter in a strong magnetic 
field, and, hence, will influence the cooling history of a 
neutron star (Reddy et al. 1999). 

We have also shown that above certain density the 
entropy for various branches of spin polarization in neu- 
tron matter with the Skyrme interaction in a strong 
magnetic field demonstrates the unusual behavior be- 
ing larger than that of the nonpolarized state. To 
clarify this point, we have provided the correspond- 
ing low-temperature analysis. It has been shown that 
this unexpected behavior should be addressed to the 
dependence of the entropy of a spin polarized state on 
the effective masses of spin-up and spin-down neutrons 
and to a certain inequality constraint on them which is 
violated in the respective density range. 

It is worthy to note that in the given research a neu- 
tron star matter was approximated by pure neutron 
matter. This should be considered as a first step to- 
wards a more realistic description of neutron stars tak- 
ing into account a finite fraction of protons with the 
charge neutrality and beta equilibrium conditions. In 
particular, some admixture of protons can affect the on- 
set densities of enhanced polarization in a neutron star 
matter with the Skyrme interaction (Isayev & Yang 
2004a). Nevertheless, at such strong magnetic fields, 
one can expect that the proton fraction is relatively 
small and even can completely disappear in the dense 
interior of a magnetar (Mao et al. 2003). 
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